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ABSTRACT 

We present high angular resolution continuum observations of the high-mass protostar NGC 7538S 
with BIMA and CARMA at 3 and 1.4 mm, VLA observations at 1.3, 2, 3.5 and 6 cm, and archive 
IRAC observations from the Spitzer Space Observatory, which detect the star at 4.5, 5.8, and 8 /im. 
The star looks rather unremarkable in the mid-IR. The excellent positional agreement of the IRAC 
source with the VLA free-free emission, the OH, CH3OH, H2O masers, and the dust continuum 
confirms that this is the most luminous object in the NGC 7538 S core. The continuum emission at 
millimeter wavelengths is dominated by dust emission from the dense cold cloud core surrounding the 
protostar. Including all array configurations, the emission is dominated by an elliptical source with 
a size of ~ 8" x 3". If we filter out the extended emission we find three compact mm-sources inside 
the elliptical core. The strongest one, Sa, coincides with the VLA/IRAC source and resolves into a 
double source at 1.4 mm, where we have sub-arcsecond resolution. The measured spectral index, a, 
between 3 and 1.4 mm is ~ 2.3, and steeper at longer wavelengths, suggesting a low dust emissivity 
or that the dust is optically thick. We argue that the dust in these accretion disks is optically thick 
and estimate a mass of an accretion disk or infalling envelope surrounding S,4 to be ~ 60M Q 
Subject headings: ISM: clouds - (stars:) circumstellar matter - stars: formation - stars: pre-main 
sequence - submillimetcr 



1. INTRODUCTION 

The giant molecular cloud southeast of the H II 
region NGC 7538 at a distance of 2.65 ± 0.12 kpc 
(Moscadclli ct al. 2009) , is a well known site of high -mass 
star formatio n (iWerner et al|1979t l Qiu et a,l.l l201lh. The 
early work by W erner et al.1 (|1979| ) showed that star for- 
mation occurs in three activity centers, each of which 
maybe in a different evolutionary stage. The IRS 1-3 
region dominates the luminosity (Lbol ~ 2 10° Lq) and 
is the closest to the H II region. Each of the three in- 
frared sources coincides with ultracompact H II regions. 
NGC 7538 harbors several young massive stars, of which 
at least three (IRS 1, IRS 9 and NGC 7538 S) drive pow- 
erful molecular outflows and appear to be surrounded by 
accretion disks. 

2 Predoctoral Fellow, Smithsonian Astrophysical Observatory 



The massive, ^25 M Q , young star IRS 1 appears to be 
the most massive m ember of a young embedded clus- 
ter (IQiu et al.ll20ll. It powers a collimated thermal 
jet dSandell et al.l 20091) and drives an ext ended molecu- 
lar outflow (ICorderl 120081: lQiu et al.ll2011IL IRS 1 is still 
heavily accreting (jCorderl 120081: IQiu et all 1201 ID . i.e., it 
must be surrounded by an accretion disk. However, no 
one has yet managed to directly detect and image the 
disk. 

NGC 7538 S, - 85" south of IRS 1 is much colder, 
has lower luminosity (Lbol ~ 1-3 10 4 Lq) and coincides 
with maser emission from OH, CH3OH, and H2O, but 
not any near-IR source, while the third region, IRS 9 
(Lbol ~ 2 10 4 L Q ), was found to be associated with a 
bright near- and mid-infrared source. Submillimeter con- 
tinuum imaging at high angular resolution starts to re- 
solve the molecular cloud into a lumpy filamentary struc- 
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ture c onnecting the individu al cores ()Sandell fc Sieversl 
[2001 IReid fc Wilsonl [200l . IRS 9 illuminates a 
prominent reflection nebulosity, has only weak free- 
free emission and drives an extreme high-velocity out - 
flow seen almost pole on (ISandell. Goss fc Wrightll2005h . 
iSandell. Goss & : Wrightl (|2005l ) found several outflows in 
the vicinity of IRS 9, suggesting that there is also a young 
cluster surrounding IRS 9. 

Figure [T] shows a large IRAC color image of NGC 7538 
with the three major star forming regions labeled as 
IRS1, IRS 9 and S, respectively. NGC7538S is a faint, 
barely visible source at the outskirt of a small cluster of 
IR sources. NGC 7538 S is not seen at 3.6 /Ltm and is not 
even particularly red due to the strong absorption from 
the surrounding cloud (Section |4~T|) . 

In this paper we focus on NGC 7538 S, which is the 
youngest high-mass formin g core in NGC 7538 an d lo- 
cated - 80" south of IRSl (|Sandell fc Sieversl 120041 ). In 
this core there is a massive y oung star, which coincides 
with an OH 1665 MHz maser. Sandcll. Wright & Forstcr 
( 2003) observed the core with BIMA in 4 mm continuum 
and in HCO+ and H 13 CN J = 1 -> with 3'.'7-resolution 
and found that it contained a high mass protostar sur- 
rounded by a large rotating circumstellar accretion disk 
driving a ver y compact molecular outfl ow. Follow-up ob- 
servations by ISandell fc Wrightl (|2010T ) confirm that the 
star is surrounded by a rotating accretion disk, smaller 
than originally thought, and that it may have fragmented 
into several protostars. 

We present new BIMA and CARMA observations from 
72 GHz to 220 GHz of the NGC 7538 S region with spa- 
tial resolutions up to ~1". These observations are com- 
bined with observations at 4.9, 8.5, 15 and 22.5 GHz from 
the VLA, and archive IRAC images (Mid-IR) from the 
Spitzer Space Telescope. 

2. OBSERVATIONS 

2.1. BIMA Observations 

Observations were obtained with the 9-antenna BIMA 
array between 2002 and 2004 in the B, C, and D ar- 
ray configurations at 72 - 75, 86 - 89, 110, and 217 - 
220 GHz. A summary of all the observations is given in 
Table [TJ which also lists the angular resolution and the 
rms-noise level achieved at each frequency. Most of these 
observations were designed to study the chemistry and 
kinematics of the accretion disk, its outflow and the sur- 
rounding cold, dense cloud core. However, in each case 
we also included one or several 100 MHz wide bands se- 
lected to be free of line emission, thereby enabling us to 
measure the true continuum emission from the disk and 
the surrounding cloud. Table Q] therefore also lists the 
main molecules and isotopes targeted for each frequency 
setting. 

At 72 - 75 GHz (A 4.1 mm) we combined the line-free 
data into a multifrequency synthesis (MFS) image. This 
image has less sensitivity than the rest of our continuum 
images, partly because not all receivers could be tuned to 
this frequency. Because of the low frequency the synthe- 
sized beam is also larger than at other frequencies (c.f. 
TableQ]). At 89 GHz (A 3.4 mm), we combined additional 
data obtained in the B, C a nd D array configurations 
with the data reported by ()Sandell. Wright fc Forsterl 
2003) to produce deep high fidelity continuum images 



and spectral line images. These observations use the 
same correlator setup and were r educed as described by 
ISandell. Wright fc Forsterl ((2001 . 

The 110 GHz (A 2.7 mm) observations were the only 
ones completely dedicated to continuum observations. 
Here we made an MFS image from observations in the B 
& C array configurations using 16 100-MHz bands span- 
ning over 3 GHz in the upper and lower sidebands of 
the 1st LO. At 220 GHz (A 1.4 mm) the only line-free 
bands were three 50 MHz windows in the lower sideband, 
but we had more tracks than at 110 GHz, and we were 
therefore still able to reach an acceptable noise level. 

The data were reduce d and imaged in a standard wa y 
using MIRIAD software (|Sault. Teuben fc Wrighil995l) . 
The quasar 0102+584 was used as phase calibrator and 
Mars and 3C 84 for flux and passband calibration. The 
uncertainty in the absolute amplitude scale is ~ 20%, 
but we were able to improve the calibration by apply- 
ing coherence corrections, see Appendix I. As a check on 
the flux density scale we also imaged IRS 1 and IRS 2, 
which are at the edge of the ~ 2 arcmin primary beam 
in the 3 mm band. In this paper we discuss only the 
continuum data, which allow us to derive an estimate of 
the mass and the physical properties of the rotating disk 
surrounding the high-mass protostar NGC 7538 S. The 
spectral line observations and the physical properties of 
the HII regions associated with IRS1, IRS 2 and IRS 3 
will be discussed in Paper II. 

2.2. CARMA Observations 

Observations were obtained with the 15-antenna 
CARMA telescope in 2009 July - August in the C, D, 
and E antenna configurations at 220 GHz, and in the B- 
configuration in 2009 December and 2010 January at 89, 
110, and 220 GHz. The CARMA observations used two 
500 MHz bands in upper and lower sidebands of LOl for 
a total continuum bandwidth of 2 GHz. The data were 
reduced and imaged in a standard way using MIRIAD 
software. 

In these observations we centered the antenna point- 
ing on IRS1 and NGC 7538 S to minimize the primary 
beam errors. Previous CARMA mosaic observations of 
NGC 7538 were designed to map the large scale molec- 
ular outflows; pointings were not centered on IRS1 or 
NGC 7538 S, and uncertainties in the primary beam pat- 
terns have prevented accurate deconvolution of the side- 
lobes of IRS1 from NGC 7538 S. Optical offset pointing 
was use d to ensure ac curate pointing of the CARMA an- 
tennas ()Corder|[20T0l) . Rapid switching between IRS1 
and NGC 7538 S was used to calibrate NGC 7538 S. The 
quasar 0102+584 was used as a gain and phase calibra- 
tor, and Uranus and MWC 349 for flux and passband 
calibration. The uncertainty in the absolute amplitude 
scale is ~ 20%. The strong compact emission from IRS 1 
was used to self calibrate the IRS 1 pointing. 

Images were made at 89, 110 and 220 GHz bands by 
combining the two 500 MHz in upper and lower sidebands 
in a multifrequency synthesis mosiac. The mosaic images 
combine data with the three different primary beams for 
the 6.1 and 10.4 m CARMA antennas. The weighted 
mean observing frequencies from the MFS synthesis are 
87.8, 111.1 and 222.2 GHz. 

2.3. VLA Observations 
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Fig. 1. — This false color Spitzer IRAC image shows the location of IRS 1, IRS 9 and S relative to the NGC 7538 H II region. The 4.5 
/im filter is coded blue, the 5.8 /^m filter is coded green and the 8.0 fim filter is coded red. The boundaries of the H II region stand out 
due to strong PAH emission from the surrounding molecular cloud. IRS 1, IRS 9 and IRS 4~ 80" northwest of IRS 1 are all saturated in 
the IRAC images. 

Our observations with the Very Large Array (VLA) of 
the National Radio Astronomy Observatory (NRAOfl 
at 3.6 cm and reduction of un published 6 cm archive 
data w ere already described in ISandell. Goss fc Wrightl 
(2005). Here we therefore only summarize the obser- 
vations and the essential parameters needed to under- 
stand and interpret the results. The 3.6 cm observations 
were done in the BnA configuration on 2003 October 
14. The on-source integration time was 1.3 hrs and the 
synthesized beam, full width half maximum (FWHM), 
was T.'22 x 0"47, (p.a. = -2°). The rms-noise is - 35 
/zJy beam -1 . The peak flux density on NGC 7538 S is 

2 The NRAO is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 



2.3 mjy beam -1 . The 6 cm observations were done in 
the B array on 1989 April 28. The on-source integration 
time was 2.7 hrs and the synthesized FWHM was 1'.'35 x 
1"09 (p.a. —24°) with an rms noise ~ 100 /iJy beam -1 . 

The 1.3 cm observations were carried out in the C con- 
figuration on 2004 March 19, pointing at NGC 7538 S. 
The integration time was 2.2 hrs. These observations 
were carried out in 4-IF continuum mode, with one 
50 MHz bandwidth pair centered on 22.485 GHz and 
one 0.195 MHz bandwidth pair was centered on the 
strongest H 2 Q maser component in NGC 7538 S, see e.g. 
IReid fc Mentenl (jl"997r i for a more thorough discussion of 
this procedure. In order to determine what maser fea- 
ture was the strongest at the time of our observations, 
we had a short test run in 2004 February, which showed 
that the strongest velocity component was centered on 
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V; sr = —53.7 km s -1 . On March 19, the flux density of 
this velocity component was 126 Jy beam -1 . The ini- 
tial calibration of the data was done using 2322+509 as 
phase calibrator and 3C48 as the flux density calibra- 
tor. We then used the strong maser feature as a phase 
calibration reference, which considerably increased the 
signal-to-noise in our final image. The synthesized beam 
FWHM for these observations was 0'.'94 x 0'.'78, (p. a. = 
— 10°). The rms-noise is ~ 0.22 mJy beam -1 . The peak 
flux density in the image is 3.3 mJy beam -1 . 

We also used a deep 2 cm observation from the VLA 
historical archive (project AM0487) conducted in the A 
array configuration on June 30, 1995. The observing 
frequency was 14.94 GHz, the bandwidth 50 MHz and 
the on source integration time 4.4 hrs. 3C 84 was used for 
flux calibration and 2229+695 for phase calibration. The 
synthesized FWHM is 0'.'139 x 0'/109 and the rms noise 
~ 80 /iJy beam" 1 . Since the observations were centered 
on IRS 1, the emission from NGC 7538 S is attenuated by 
about a factor of two by the primary beam with a peak 
flux density of 0.4 mJy beam -1 . 

2.4. IRAC archive data from Spitzer 

We have retrieved and analyzed images of NGC 7538 
taken with the Infrared Array Camera (IRAC). IRAC is 
the mid-infrared camera on the Spitzer Space Telescope 
with four arrays at 3.6, 4.5, 5.8, and 8.0 /im. Each ar- 
ray has 2 56 x 256 pixels with the same plate scale, 1'.'22 
per pixel ()Fazio fc M egcath 2004) . The observations dis- 
cussed here were obtained on December 23, 2003 as part 
of the GTO program P201 (G. Fazio) and observed in the 
High Dynamic Range (HDR) mode with a three point 
dither for each pointing. In the HDR mode images are 
taken with a short (0.6 s) and a long (12 s) integration 
to capture b oth bright and fain t sources in each point- 
ing. IRS 11 (|Werner et al.| [T979) is the brightest source 
in the NGC 7538 S field and it is saturated in the long 
integration exposures at 8 /im, resulting in bandwidth 
effect trailfl At 8 /im we therefore only use the short 
exposure images. 

We processed both the short (0.6 s) and the long (12 s) 
integration basic calibrated data (BCD) IRAC frames in 
each channel using the artifact mitigation software devel- 
oped by Sean Carey and created mosaics using MOPEX. 
The saturation of the 8 /im images is particularly severe 
and could not be corrected for by the available tools. At 
8 /im we therefore used only the short integration im- 
ages to identify sources and derive photometry. Because 
the whole NGC 7538 complex is associated with exten- 
sive nebulosity and strong PDR emission (Figure Q}, all 
photometry was derived using PSF fitting with APEX, 
which is part of the MOPEX reduction package. For 
sources which APEX failed to detect at one or several 
wavelengths, we used the APEX user list option to sup- 
ply the coordinates for the source to successfully derive 
a PSF fit. 

We also retrieved and analyzed Spitzer MIPS images 
of NGC 7538. NGC 7538 S is not saturated at 24 micron, 
but IRS11 is and NGC 7538 S is blended with the satu- 
rated PSF from IRS11, so it is not possible to derive a 

3 see the Infrared Array Camera Data Handbook , 
|http://ssc. spitzer. caltech.edu/irac/ /dh/irac datahandbook3.0.pdf 
lor details. 
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Fig. 2. — High-resolution images of NGC 7538 S made with 
CARMA in B array configuration at 222 GHz (gray scale), 111 
GHz (cyan contours), and 87 GHz (green contours), superimposed 
on the deep image of the NGC 7538 S cloud core with CARMA at 
224 GHz in the R, C, D and E array configurations (red contours). 
The beam size (FWHM) for the B-array images at 222, 111 and 
87 GHz are shown at bottom left. Contour intervals: 13.3 mJy at 
222 GHz, 3, 6, 9 mJy at 110, and 2, 3, 6, 9 mJy at 87 GHz 

flux. At 70 micron the whole region from IRS1 south- 
ward to IRS11 and NGC 7538 S is saturated and contains 
no information. 

3. RESULTS 
3.1. CARMA 

Figure [5] shows the deep 222 GHz B, C, D, and E array 
image of the NGC 7538 S cloud. The emission is domi- 
nated by the elliptical core, which at the western side 
has an extension to the northwest and a curved filamen- 
tary like structure to the southwest. The northwestern 
extension has one embedded sub- mm source, the H2O 
maser IRS 11S, which we detected at 111 and 222 GHz, 
but not at 88 GHz. In Figure [2] this source is seen as 
only one cyan contour level. The sa me source was also 
marginally de tected bv lCordeil ((2008), who labeled it S/. 
ICorderl (|2008f) also found two embedded protostars in the 
fainter core in the southwestern filament, which is labeled 
as BIMA source South in Table 2. The core is clearly seen 
in the deep 222 GHz image, but we did not have enough 
sensitivity in the B-array images to recover any embed- 
ded sources. At the highest resolution, 0'.'2, i.e. the 222.2 
GHz B-array image (Figure [5] gray scale and Figure 8), 
the high-mass protostar, S^, has been resolved into two 
components: Sai> and S^2- The flux densities of the 
components were estimated from Gaussian fits. The po- 
sitions, integrated fluxes, deconvolved sizes, and position 
angles are given in Table 2. 

3.2. VLA 

The VLA observations show a double source at 4.9, 
8.5, and 22.5 GHz (Figure [3]). The stronger component 
(a) is centered on NGC 7538 S, while the second compo- 
nent (b) is ~ 1'.'5 to the north west at a position angle of 
—34°. Both components appear elliptical, and Gaussian 
fits, when the components are sufficiently resolved (Table 
2), show that their major axis is aligned in the same di- 
rection, suggesting that the two components are part of a 
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Fig. 3. — VLA continuum images centered on NGC 7538 S at 
1.3, 3.6, 2 and 6 cm. The 2 cm image (grey-scale only) had IRS 1 
as pointing center and has poor SNR because it is in the outskirts 
of the primary beam and has very high angular resolution, ~ O'.'l. 
The integrated intensity, however, agrees with the flux densities 
measured at other frequencies and confirms that the free-free emis- 
sion originates in a highly collimated ionized jet. Position offsets 
are w.r.t. NGC 7538 S a 

collimated jet emanating from NGC 7538 S into the blue- 
shifted molecular outflow lobe p owered by NGC 7538 S 
(|Sandell. Wright k, Forster|[2b0 3V The high angular res- 
olution image at 14.9 GHz confirms that the main com- 
ponent has the appearance of a highly collimated jet, 
although it does not have enough sensitivity to clearly 
detect the northwestern component, which in this picture 
would be the tip of the jet. Our deepest image, at 8.5 
GHz, reveals a faint extension to the southwest, which we 
interpret as the counter jet traveling into the red-shifted 
outflow. A least squares fit to the sum of both compo- 
nents yields a spectral index of 0.14, which confirms that 
we are looking at win d-ionized free-free emission from a 
collimated jet, see e.g. I Reynolds! (|1986D . 

3.3. BIMA 

Images of NGC 7538 S at 74, 89, 109, and 217 GHz are 
shown in Figure 2] The 109 and 217 GHz images appear 
very similar and agree quite well with the deep CARMA 
image at 222 GHz, although the CARMA image recovers 
more of the extended emission in the cloud core (Fig. [2]). 
Both the 109 and the 217 GHz image show extended 
dust emission with a size of ~ 20". The emission is dom- 
inated by an elliptical core centered ~ 1" - 2" southwest 
of NGC 7538 S with a position angle of - 50°. They 
also show two secondary peaks, one ~ 5" west and 1" 
north, the other - 5" west and 6" south of NGC 7538 S. 
The former coincides within errors with an H 2 maser, 
which is also detected in all IRAC bands from 3.6-8 
fim, see Section |4~T1 The H2O maser was also detected 
as source D in the CARMA B-array observations (Table 
2). T he southweste rn core has no maser or IR counter- 
part. iCorderl (2008) identified two 3 mm sources in this 
core, but we have not been able to confirm them. The 
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Fig. 4.— BIMA continuum images of NGC 7538 S at 1.4, 2.8, 3.4, 
and 4.1 mm in grayscale overlaid with contours. The red contours 
on each image shows the free-free emission at 3.6 cm. On the 1.4 
mm image we have additionally marked the position of the two 
IRAC mid-IR sources coinciding with NGC 7538 S and the western 
H2O maser source IRS 11 S (Table 3). For each image we have 
plotted 6 linearly spaced contours between ~ 3-cr rms level and 
peak flux density. The beam FWHM is indicated at the lower left 
of each panel. Position offsets are w.r.t. NGC 7538 Sa 

image at 74 GHz looks different. Here the emission is 
much fainter than at the higher frequencies and the spa- 
tial resolution is poorer. The difference in morphology is 
most likely due to poor image fidelity and sensitivity. 

4. DISCUSSION 

4.1. The embedded stellar population - IRAC mid-IR 

imaging 

Table 3 lists the positions of the 14 sources which were 
detected at 8 /mi within a 25" radius around NGC 7538 S. 
NGC 7538 S is detected at 4.5, 5.8, and 8 /xm, but not at 
3.6 /im. There is a small cluster of four sources centered 
on IRS 11, and another faint source east of NGC 7538 S. 
One of these sources, ~ 5" west of NGC 7538 S coincides 
with an H2O maser and is also seen as a compact dust 
emission source in the BIMA continuum images at 109 
and 217 GHz. 

In Figure 6 we plot all the 14 sources from Table 3 in an 
IRAC color-color plot. Six of the 14 stars in Table 3 have 
low color indexes ([5.8]-[8.0] < 0.4 and [3.6]- [4.5] < 1.5) 
and are most likely reddened field stars. Two of them, 
J23134708+6127090 and J23134160+6126427, however, 
have [5.8]-[8.0] colors close to 0.4 and could therefore be 
more evolved pre-main sequence stars. The group of five, 
possibly six sources if we include J23134708+6127090, 
which lie to the northeast of NGC 7538 S all have [5.8]- 
[8.0] colors > 0.7, and are clearly young pre-main se- 
quence stars. The brightest star in the group is IRS 11, 
which does not have much foreground extinction, [3.6]- 
[4.5] ~ 0, but IRS 11 has a strong 8 /itm excess, and is 
clearly a young star. IRS 11 S, the H2O maser, appears 
to be much more deeply embedded than IRS 11, and it 
has a similar [5.8]-[8.0] color as IRS 11. Just north west 
of this group is J23134459+6127148, which is surrounded 
by a prominent nebulosity at 5.8 and 8 /im and is almost 
certainly a young object, but it may be associated with 
IRS 1 rather than with the NGC 7538 S cloud core. With 
NGC7538 S we therefore see a young cluster with at least 
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Fig. 7. — Integrated flux densities of NGC 7538 S as a function 
of frequency. These have been corrected for coherence loss (see 
text) and corrected for free-free emission. The solid lines show the 
spectral index, a, determined from the BIMA and CARMA data. 
The spectral index appears to steepen at longer wavelengths. See 
text. 



Fig. 5.— The detailed image of NGC 7538 S region in the 3.6, 4.5, 
5.8 and 8.0 fim. The identified point-like sources are indicated with 
crosses and circles. The isolated sources labeled by crosses have 
relatively low color index and are likely to be field star or evolved 
class III stars. The sources plotted as open circles have large color 
indeces, and are all young pre-main sequence objects. NGC 7538 S 
and IRS 11 S, both of which were detected in our CARMA images, 
are very deeply embedded ([3.6]-[4.5] > 2). NGC 7538 S was not 
even detected at 3.6 fim. 
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Fig. 6.— The color-color diagram of [3.6]-[4.5] vs. [5.8]-[8.0] for 
the 8 fim sources in NGC 7538 S. The indications of the crosses and 
circles in this figuare are the same as Figure 5. 

six sources. Our deep 1 mm CARMA images show at 
least two additional sub-mm sources, which must be in 
an even earlier evolutionary stage, or more deeply cm- 
bedded, since they have no mid-IR counterpart. 

4.2. Millimeter wavelength dust emission 

The VLA observations fSection l3.2l) show that the free- 
free emission in the 3 mm band is predicted to be ~ 
8.5 mjy and therefore the emission at mm-wavelengths 
is completely dominated by dust emission. The free- 
free corrected flux densities obtained with BIMA and 
CARMA for the elliptical core are plotted in Figure [7] as 
a function of frequency. A least squares fit to the BIMA 
89, 109 and 217 GHz data points give a spectral index, 
a ~ 2.45, which is surprisingly low. At 89 GHz, however, 



the elliptical core appears more centrally condensed and 
smaller (Table 2, Figure 0]), i.e. the observations are 
not sensitive enough to recover the cold outer parts of 
the core. Since the 89 GHz observations do not sample 
the same volume of dust, it should therefore be excluded 
from the fit. Using only the flux densities from the best 
images, at 109 and 217 GHz, gives a = 2.3 ± 0.2. This 
corresponds to an apparent dust emissivity, /3 ~ 0.4, if we 
assume that the dust is isothermal with a dust tempera- 
ture of 40 K. This is much lower th an the dust emissivity 
derived by Sandcll fc Sieversl (|2004[) from multi frequency 
single dish observations and appears suspiciously low. 

We th erefore also ex amined the CARMA data pre- 
sented in lCordeil (|2008l ) as well as the new CARMA data 
discussed in this paper. Corder finds a bigger size for the 
elliptical core, ~ 12" x 7", and we therefore cannot di- 
rectly compare his results to our BIMA results. Instead 
we compare the quoted total fluxes for the NGC 7538 S 
cloud core, which is 0.218 ± 0.007 Jy at 91.4GHz and 
0.407 ± 0.016 Jy at 108.1GHz. For the deep BCDE- 
array CARMA image at 222.2 GHz we find a total flux 
for the core of 1.9 ± 0.2 Jy. A least squares fit to all 
three data points give a = 2.5 ± 0.4, which is very sim- 
ilar to what we found from the BIMA data. If we only 
use the 108 and 222 GHz we get a = 2.2 ±0.4. Both the 
CARMA and the BIMA data (Fig. 7) show a steepen- 
ing of the spectral index at long wavelengths. If we only 
use the data points at 91 and 108 GHz, we get a = 3.8 
± 0.5, which is consistent with thermal emission from 
dust. What these data therefore suggest is that the el- 
liptical core may be so extreme, that the dust is already 
becoming optically thick at 1 mm. It is possible that 
what we observe is partly an instrumental effect, i.e., we 
filter out more emission at the highest frequencies. How- 
ever, the same trend is seen in both our BIMA and our 
high fidelity CARMA data, which at 222.2 GHz have ex- 
tremely good uv-coverage. This suggests that the trend 
is real and not just an artifact from different response 
in spatial sensitivity for an aperture synthesis array at 
different frequencies. 

One would in fact expect the surrounding cloud 
core to have a steeper spectral index than the disk, 
because most studies find that the emissivity in- 
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Fig. 8. — High spatial resolution images of NGC 7538 S ob- 
tained with CARMA at 88 (black), 111 (blue) and 222 (red) GHz. 
The synthesized beam FWHM shown in the lower left corner are 
0.87" xO. 68", 0.61"x0.52", and 0.34" x0.27"at 88, 111 and 222 
GHz respectively. Contour intervals: 2.56 mjy/beam at 222 GHz, 
2, 3, 6, 9 mjy/beam at 111 and 88 GHz. 



dex ~ 1.5 - 2 (IWrightl U99l IMasi et all U995t 
iGoldsmith. Bergin fc Langerl I1997D i.e. corresponding 
to a ~ 4. The most likely explanation to the observed 
low spectral index of the surrounding cloud core is that 
the BIMA array filters out much more of extended emis- 
sion at 21 7 GHz than at 109 GHz. The single dish obser- 
vations by Sa ndell fc Sieversl (|2004[) shows that the cloud 
core surrounding NGC 7538 S has a radius of ~ 30", while 
our interferometer observations find a core radius of less 
than half of this value, suggesting that most of the cloud 
core is resolved out in the BIMA observations. Both 
the 217 and 109 GHz observations recover spatial scales 

of ~20"; therefore the results for the disk surrounding 
NGC 7538 S should not be affected. 

Although variations in the observed spectral index or 
dust emissivity are usually interpreted as changes in 
dust composition, such changes can also be caused by 
temperatur e and density gradients if the dust is opti- 
cally thick. iBeckwith fc Sargent! (|199lD found that most 
T Tauri stars have spectral indeces between ~ 2 - 3, 
which they showed to result from T Tauri stars hav- 
ing optically t hick in ner disks at millimeter wavelengths. 
iDutrev et al.l (1996) did a more careful study using the 
IRAM interferometer at 2.7 mm and confirmed that all 
the stars in her sample had optically thick inner disks. 
It is now also well established that disks around more 
evolved T Tauri and Herbig Ae stars have low dust emis- 
sivity, but for these stars the drop in emis sivity is due 
to gr ain growth, not optically thick dust (|Natta et al.l 
[200l . 

For high-mass protostellar objects, single dish surveys 
sho w a much larger scatter in the observe d dust emissiv- 
ity. iWilliams. Fuller fc Shridharanl (pOOl did a SCUBA 
850 and 450 /jm survey of 68 protostellar objects and 
found a mean spectral index, a = 2.6 ± 0.4, correspond- 
ing to a dust emissivity, /3 = 0.9 ± 0.4, suggesting that 
the dust emissivity in high-mass protostellar objects has 
largely the same characteristics as dust disks surround- 
ing T Tauri and Herbig Ae/Be stars. On the other hand, 
iHill et al.l(|2006f ). who did a similar survey using SCUBA 
observations at 850 and 450 fim and SIMBA observations 
at 1.3 mm of 212 cold massive cloud cores, argue that 
their results are in agreement with a typical /3-value of 
2, although their results show considerable scatter in the 
derived /3 indeces to both higher and lower f3- values. 
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Fig. 9. — Overlay of H2O maser positions on the VLA 3.6 cm 
continuum image. The contours for the continuum image are lin- 
ear with seven contours between 0.10 mjy beam -1 and the peak 
flux, 2.47 mjy beam -1 . The H2O maser spots are marked with 
triangles, where the red triangles show red shifted maser spots, the 
blue ones show blue shifted masers, and the black ones are near the 
sy stemic velocity of the cloud. The filled circle shows the position 
of the IRAC mid-infrared sourc e. The filled star symbols are RCP 
(right circularly polarized) OH 1665 maser s pots, while t h e ope n 
circles are LCP OH 1665 maser spots from Ergon et all (2000), 
which are measured to an accuracy of Of.' 1. The cross marks th e 
position of the CH3OH class II maser from[Pcstalozzi ctal. (2006), 
who quote an astrometric accuracy of / /5. Position offsets are 
w.r.t. NGC 7538 S a 

However, all published millimeter array observations of 
high-mass protostellar objects find consistently low /3 val- 
ues. For the accretion disk in IRAS 20126+4104, which is 
barely resolved in mm-continuum, iCesaroni et al.l (|2005[ ) 
find a spectral index of 2.7 from a least squares fit to 
all published Plateau de Bure Interferometer (PdBI) 
data at frequencies 86.6, 90.6, 96.5, 217.6 and 239.6 
GHz, the Owens Valley Radio Observatory (OVRO) 
millimeter wave array data at 113.2 and 228.1 GHz, 
and a 7 mm (42.8 GHz) da t a poi nt from the VLA. 
iBeuther. Schilke fc Wvrowskl (2004) were able to re- 
solve the high mass protostellar core mml in IRAS 
20293 + 3952 with PdBI at 98.0 and 244.9 GHz, and 
found a spectral index of ~ 3.4 at the outer edge of the 
core decreasing to values as low as ~ 2.1 at the very 
center of the core. SMA observations of the massive pro- 
tostellar core I RAS 18089 - 1732 w ith the SMA at 217 
and 354 GHz (jBeuther et al.l 120041 ) give similar results 
from analysis of fluxes in the u-v plane; j3 ~ 1.2 for large 
spatial scales, 5" - 10", and decreasing to ~ 0.5 for small 
spatial scales, 2'.'7 - 3"6. Therefore, in the few cases 
where high-mass protostellar objects have been resolved, 
they all indicate that the dust emissivity drops towards 
the center of the core, which is believed to be dominated 
by emission from an accretion disk. 

The high-mass protostellar objects discussed here are 
very young. It is therefore unlikely that grain growth 
would have had time to occur. A more plausible explana- 
tion for our low spectral index between 111 and 222 GHz 
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is that the disks are optically thick at mm-wavelengths. 
Our results on NGC 7538 Sa strongly favor high optical 
depth, rather than a population of large grains in the 
disk. It is clear from the BIMA images (Figure |4j that 
Sa starts to dominate the emission of the elliptical core 
at the longest wavelengths, suggesting that one can now 
see the center of the accretion disk, which is expected 
to be hotter and denser than the surrounding infalling 
envelope. The center of the disk therefore dominates the 
emission. 

If the disk is optically thick at 111 and 222 GHz, as 
our observations suggest, one would really need to model 
the dust emission in order to derive a disk mass. How- 
ever, we can get an estimate of the disk mass by assum- 
ing that the dust emission is optically thin and roughly 
isothermal at the longest wavelength for which we have 
a good high resolution image, i.e. at 89.3 GHz. The 
total mass of gas and dust, M, can then be expressed 
as M = S l/ D 2 _/(K l/ B l/ (T d )), where B„(T d ) is the Planck 
function, T d is the dust temperature, D is the distance, 
and S„ is the integrated flux density at the frequency 
v. Because of the unce rtainty in dust em issivity, we 
do not want to use the iHildebrandl (|1983t l dust opac- 
ity and extrapolate it from 250 /im. Instead we have 
ch osen to use the dust mass o pacity, k as tabulated 
by lOssenkopf fe Hennin g (1994) at 1.3 mm. Their re- 
sults suggest a dust mass opacity Ki.3 mm — 0.01 cm 2 g _1 
for gas densities of 10 s cm -3 and grains with thick ice 
mantles, assuming the ratio of gas to dust mass is 100. 
These are the type of dust grains and densities that 
one would expect to have in the NGC 7538 S accretion 
disk. At frequency u, the dust mass opacity is there- 
fore k v = 0.01 x (230/i/)-" cxd?g~ x . At 91.4 GHz we 
now obtain a dust mass opacity of ~ 0.0025 cm 2 g _1 
if we assume ft = 1.5, which is probably appropriate 
for this case. If we further assume a dust tempera- 
ture of 35 K, see below, a distance 2.65 kpc, and cor- 
rect the observed flux density at 91.4 GHz for free- free 
emission (~ 8.4 mJy) we obtain a disk mass ~ 60 M . 
T his is about two times small e r tha n what was reported 
in iSandell. Wright fc Forsterl J2003), but in reasonable 
agreement with th e dynamical mass derived from H 13 CN 
J=l->0, 20 M Q (jSandell. Goss fc Wrightl[2005h . Since 
the size of the rotationally supported part of the disk is 
smaller than what we measure in continuum, we would 
expect the mass derived from continuum to be higher. 

The assumption of a dust temperature o f 35 K ap- 
pears .justified, see e.g. ISandell fc Sieyer jjpOOl and 
ISandell fc Wright] (pOlOh . IZheng et all (1200 lh find am- 
monia to be optically thick in the direction of NGC 7538 S 
with a temperature of 25 K. The temperature of the 
disk is likely to be less than 50 K, because then 
we would exceed the obs erved far infrared luminosity 
|Th ronson fc Harperll 19791) . The dust temperature must 
be higher than 25 K, the temperature of the surround- 
ing cloud core, since ISandell fc Wright! pOlCh see some 
self-absorption in virtually every molecule seen in the 
disk. We can get another upper limit by looking at the 
analys is of the CH 3 CN transitions by ISandell fc W right 
(2010), where they derive a rotation temperature = 52 
± 10 K. Since the low K-transitions are optically thick 
(jSandell fc Wright 2010j), this leads to an overestimate of 
the rotational temperature. Contribution from the hot 



outflow will also raise the temperature. Therefore a dust 
temperature of 50 K appears to be a reasonable upper 
limit. If we were to increase the dust temperature to 50 
K, this would change our mass estimate by less than a 
factor of two. 

We can also do an order of magnitude estimate of the 
disk mass if we assume that the inner disk is optically 
thick at 1.4 mm using a geometrically thin disk model, 
which has been successfully used to an alyze millimeter 
continuum emission from T Tauri stars feeckwit h et al.l 
119901 : [Dutrev e t al. 1996J). For an optically thin disk the 
optical depth, r„ = k v x E(r) / cos(9) , where k„ is dust 
mass opacity, S(r) is the surface density as a function of 
radius, r, and 9 is the inclination angle of the disk (9 — 
90° for an edge on disk). The surface density is assumed 

T _ 

to be a power law of radius, i.e. S(r) = S D x ( — ) p . 

Ro 

If we define the disk radius, Rd, as the radius where the 
optical depth, t v = 1, we can express the disk mass Md 
as 



M D = -^^i? 2 

2 - p K„ 



(1) 



which, if we adopt a dust mass opacity of 0.009 cm 2 g 1 
at 217.3 GHz gives: 



Md 



78.6 



P 



cos{9)( 



loooAir 



[M G 



(2) 



For an optically thick radius of 0'.'5 (1400 AU), roughly 
the size measured with CARMA, a relatively flat mass 
distribution, i.e, p = 1, and a disk inclination of 70° , 
we derive a mass of ~ 50 M Q for the optically thick 
disk component, which is similar to what we estimate 
assuming the dust is optically thin at 91.4 GHz. The 
radius of the optically thick disk is probably somewhat 
smaller, so this is in fact an upper limit. In order to be 
able to determine the radius of the optically thick inner 
disk, it would be necessary to have observations with 
a resolution of O'.'l - 0'.'2 at at least three frequencies, 
one of which would need to be below 100 GHz. Such 
observations are now completely feasible. 

We can estimate the dust opacity towards S^ from 
our high resolution CARMA data. We need to subtract 
the emission from the thermal jet, which from the fit to 
our VLA data is predicted to be: 8.3, 8.6, and 9.5 mJy 
at 87.9, 111.1, and 222.2 GHz respectively. Fitting the 
three data points at the above frequencies on S^ after 
subtracting the contribution from the free-free emission 
(here we add the flux for S^i and S^2 at 222.2 GHz 
,and the flux is assumed to be 47 mJy), we obtained 
a spectral index a — 2.1 ± 0.3, assuming that all data 
points have a 10% error. This is a strong indication that 
the dust emission towards Sa, our high mass protostar, 
must be optically thick. Because the dust emission from 
the NGC 7538Sa disk is rather strong, it should now be 
possible to detect the dust emission with the eVLA both 
at 22 and 43 GHz, which would allow us to get a better 
estimate of the dust emissivity index. 

4.3. Protostars NGC 7538 Sa, S b , and S c 

In addition to NGC 7538 Sa, the high-resolution 
CARMA image (Figures [2] and [8]) reveals two compact 
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components (Sb and Sc) southwest of Sa along the ridge 
of the accretion disk. No radio continuum emission from 
either of these two components has been detected with 
the VLA, suggesting that they are in an early protostellar 
phase. They may have been produced from fragmenta- 
tion due to an instability of the material in the accretion 
disk of NGC 7538 S. They are the youngest members in 
the young stellar object (YSO) group of NGC 7538S in- 
cluding IRS sources detected with Spitzer. When the 
elliptical core started collapsing, Sa was probably the 
first protostar to be surrounded by a rotating disk which 
already at that point was very dense and optically thick. 
We know that there is a velocity gradient over the whole 
8" elliptical core. Whether Sb and Sc fragmented from 
the massive, unstable disk around Sa or whether the core 
fragmented into three protostars is unclear. It would be 
important to know but it is beyond the scope of our cur- 
rent paper. We see a very dense rotating core in which 
we find three extremely young protostars. Sa has the lu- 
minosity and mass to form a star which may evolve into 
an early B or late O star. We do not know the luminosity 
of the other two, but judged on the Amm flux (and they 
have no free- free), they are similar or slightly lower mass, 
so they probably will form mid- to late B-stars. The 
Bonner- Ebert mass estimated from the observed size and 
flux density at 225 GHz for dust temperatures 25 to 50 
K are 0.01 to 0.1 M Q . The mass of the compact sources 
exceeds the Bonnor-Eb ert mass suggesting that t hey are 
local accretion centers ()Bonnorl ll956; Ebern ll957| ) . Com- 
petitive accretion will make Sa grow faster, so it will al- 
most certainly be the most massive star. The accretion 
flow, withi n the errors, is centered on what we now know 
to be S A (ISandell fc Wrightl lUlCft . As we have shown 
in this paper this is also the most luminous and massive 
protostar in the elliptical core. It is therefore reasonable 
to assume that accretion will continue until Sa forms an 
HII region, which then will halt the accretion onto the 
star. However, when this happens, and after the HII 
region expands, it will also halt the accretion onto the 
nearby protostars, so Sa will remain the most massive 
star. 



5. CONCLUSION 

We used BIMA and CARMA observations of NGC 
7538S at millimeter wavelengths, VLA observations at 
centimeter wavelengths, and Spizter in the mid IR bands, 
to study the high-mass star formation in this YSO group. 



The observations at millimeter wavelengths show that 
the continuum emission from NGC 7538 S is dominated 
by dust emission from an 8" x 3"source, elongated north- 
east to southwest, suggesting an accretion disk. The ob- 
served spectral index of a ~ 2.3 implies that the dust in 
the disk either has a low dust emissivity or is optically 
thick. 

The elliptical core breaks up into three Amm sources. 
The brightest one, NGC 7538 Sa, is extended and is pos- 
sibly a binary. This Amm source coincides with an OH 
1665 MHz maser, a Class II methanol maser, a cluster of 
H2O maser spots, a faint free-free VLA jet, and a heav- 
ily obscured mid-IR source (Spitzer IRAC and IRS data). 
It drives a very compact bipolar outflow and appears to 
be surrounded by a rot ating disk with a radius of ^2" 
(jSandell fc Wrightl |2010|). The other two A3 mm sources 
along the central ridge of the clump both appear point- 
like, alt hough one of them may drive another compact 
outflow rtCorderll2008ft . 

NGC 7538Sa appears to be the most powerful source 
in this group and produces an ionized bipolar outflow, 
detected by the VLA, in a direction perpendicular to the 
major axis of the disk. At millimeter wavelengths, with 
sub-arcsec resulution, CARMA has resolved this source 
(Sa) into two sub-components distributed along the ma- 
jor axis of the accretion disk. Located southwest of Sa 
along the major axis of the accretion disk, two compact 
millimeter sources Sb and Sc are also detected at 222 
GHz with CARMA, suggesting that they may be the re- 
sults from the fragmentation of the accretion material 
in NGC 7538S. Along with the highly reddened objects 
detected with Spitzer at mid IR bands, the protostars 
detected at millimeter wavelengths suggest that NGC 
7538S consists of a group of very young stellar objects 
in the NGC 7538 complex. 

The BIMA array was operated by the Universities of 
California (Berkeley), Illinois, and Maryland with sup- 
port from the National Science Foundation. Support 
for CARMA construction was derived from the states 
of California, Illinois, and Maryland, the Gordon and 
Betty Moore Foundation, the Kenneth T. and Eileen L. 
Norris Foundation, the Associates of the California In- 
stitute of Technology, and the National Science Foun- 
dation. Ongoing CARMA development and operations 
are supported by the National Science Foundation under 
a cooperative agreement, and by the CARMA partner 
universities. 



APPENDIX 
ATMOSPHERIC DECORRELATION 

The resolution and dynamic range of millimeter wavelength aperture synthesis images are limited by atmospheric 
turbulence. Atmospheric phase noise reduces the observed amplitudes by decorrelating the signals received at each 
antenna. In an analysis of data obtained with the Hat Creek interferometer at 86 GHz we found a baseline dependent 
phase structure function with a power law index, f3 between 0.6 and 1.7. A slope close to the Kolmogorov 2D value 
0.67 is obtained for longer baselines and stable weather conditions; the slope approaches the Kolmogorov 3D value 
1.67 on short baselines and in turbulent weather. The rms pathle ngth on a 1 km baseline ~ 1 mm, with a variation 
by a factor 4 over a few days even during good weather conditions iWrightl (|1996h . 

We estimated the atmospheric decorrelation for the NGC 7538 observations by analysing the phase calibrator data. 
The phase calibrator was observed at ~ 25 minute intervals and samples the same atmospheric fluctuations as the 
observations of NGC7538. Although the short calibrator observations do not have sufficient signal-to- noise to de- 
termine the radio seeing for each observation of NGC7538, we can estimate the seeing over several hours from the 
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observed calibrator phase scatter. Atmospheric phase noise increases the measured FWHM of the calibrator image. 
We calculated amplitude corrections for atmospheric decorrelation by fitting a circular Gaussian model to images of 
the phase calibrator. 

The BIMA data for NGC 7538 were corrected by applying an amplitude correction for atmospheric coherence, by 
specifying the FWHM in arcsec, of a circular Gaussian fit to the phase calibrator(s) for the target source. The derived 
scaling factors from decorrelation loss were 1.11 at 74 GHz, 1.08 at 89 GHz, 1.50 at 109 GHz, and 1.18 at 217 GHz. 
The flux densities in Table 1 have been corrected for these decorrelation losses. 
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TABLE 1 

BIMA OBSERVING LOG. 



Array Configuration rest frequency 
[GHz] 


Bandwidth 
[MHz] 


Synthesized beam 
" X " pa = ° 


rms 
[mjy beam -1 ] 


Observing dates 


Frequency setting: DCN J=l— >0 


BCD 73.9 


400 


9.1 x 7.2 78.9 


3.7 


22 Jul & 30 Oct 03, 07 Jan 04 


Frequency setting: HCO+ J=l->0, H 13 


CN J=l->0, 


NH 2 D In ->■ lot, SO 3=2, 2-kL, 1 




BCD 89.3 


75 


4.1 x 3.9 8.8 


3.0 


27 & 31 Jan ,& 03 May 02 

09 Aug, 04 Nov, 11, 12, & 31 Dec 03 


Frequency setting: Continuum 










BC 109.0 


800 


3.5 x 3.2 89.5 


0.54 


27 Mar & 27 Apr 04 


Frequency setting: DCN J=3->2, 13 CO J=l-rt), CH 3 CN J=12->11 






BC 217.3 


150 


2.9 x 2.4 78.4 


3.0 


13 Oct 02, 19 Apr, 1 May 03 
10, & 25 Oct 03, 3, & 5 Jan 04 
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TABLE 2 

Continuum positions and flux densities for sources in the NGC 7538 S core from VLA, 

BIMA and CARMA observations. 



Source 


Frequency 


a(J2000.0) 


<5(J2000.0) 


Ba X 6»f, a 


p. a. 






[GHz] 


m sj 


["'"] 


r"x"i 
i a j 


f°l 
I 1 








VLA Observations 








IN V..' V_ 1 OOO O cL 




23 13 44.992 


+61 26 49.6 




-f-±o± 


q fin _i_ n i & 






23 13 44.817 


+61 26 50.9 




+119 


1.76 ± 0.17 




8.46 


23 13 44.971 


+61 26 49.7 


n r v 1 


+137 


Q OK _L n HQ 


b 




23 13 44.822 


+61 26 50.8 


0.9 x 0.4 


+136 


1 f.Q _L f) 1 1 


c 




23 13 45.066 


+61 26 49.1 


i n v i 

l.U A ^U.-L 


4-117 

-i-n i 


u.yo ztz u.iz 




14.9 


23 13 44.953 


+61 26 49.8 


zL^y 94 

U.*±OA 


—36 


4.4 ± 0.4 a 


a 


22.49 


23 13 44.978 


+61 26 49.7 


0.8 x 0.1 


-37 


4.77 ± 0.06 


b 




23 13 44.797 


+61 26 51.2 


1.2 x 0.4 


-29 


1.98 ± 0.08 






BIMA 


Observations 








INLrU /Ooo b 


74.0 


23 1 3 44 836 


+61 26 50.5 


point-source b 




on Q_L/i n 




CO Q 

oy.o 


23 1 3 44 826 


+61 26 49.1 


7.4 x 3.4 


+ ( u 






i no n 


23 1 3 44 91 6 


+61 26 48.8 


8.2 x 3.9 


1 A Q 


zzUztzU 




217.2 


9^ 1 ^ Ad QDQ 


+61 26 48.8 


7.8x2.8 


+54 


1 non_L i nn 
1U2U±1UU 


South 


109.0 


9^ 1 ^ A4 90 


+61 26 42.2 


5.1x2.8 




22i4 




217.2 


23 13 44.28 


+61 26 42.2 


5.5x0.1 




96+10 


West 


i no n 

ruy.u 


23 13 44.12 


+61 26 50.2 


6.7x2.7 


-\-ZO±i- 


on 1 Q 




217.2 


23 13 44.18 


+61 26 50.0 


4.3x2.8 


+57±3 


155+20 






CARMA B- 


-array Observations 






SA2 


222.2 


23 13 45.057 


61 26 50.1 


0.49x0.33 


35 


26.1 


Sai 


222.2 


23 13 44.986 


61 26 49.8 


0.51x0.34 


-22 


40.9 


Sb 


222.2 


23 13 44.750 


61 26 48.8 






36 


Sc 


222.2 


23 13 44.510 


61 26 47.8 






34 


Sd 


222.2 


23 13 43.88 


61 26 50.2 






5 


Sa 


111.1 


23 13 45.005 


61 26 49.8 


0.61x0.58 


48.7 


20.1 


Sb 


111.1 


23 13 44.750 


61 26 48.8 






5 


Sc 


111.1 


23 13 44.510 


61 26 47.8 






8 


Sd 


111.1 


23 13 43.878 


61 26 50.2 


1.79x0.68 


20.1 


17.5 


Sa 


87.9 


23 13 45.005 


61 26 49.8 


0.84x0.37 


85.5 


13.3 


Sb 


87.9 


23 13 44.750 


61 26 48.8 






3 


Sc 


87.9 


23 13 44.510 


61 26 47.8 






2 


Sd 


87.9 


23 13 43.88 


61 26 52.2 









a At the edge of the primary beam (correction factor ~ 2). Flux density underestimated because the obser- 
vations are not sensitive enough to pick up faint extended emission in a 0'! 1 beam. 

k Compact or marginally resolved point-like source in 9" beam embedded in faint extended emission. 



High-mass protostar NGC 7538 S A 



13 



TABLE 3 

Positions and flux densities of 8 /ai IRAC sources in the vicinity of NGC 7538 S. 



Source 


a 


(2000.0) 


5(2000.0) 




S(3.6 /im) 


S(4.5 Atm) 


S(5.8 


/zm) 


S( 


1 fim) 


Comment 








h m sj 


[° ' "] 




[mjy] 


[in 


Tvl 

J yj 


[mjy] 


[mjy] 




J23134726+6126339 


OQ 


1 Q A *7 OQ 
lo 4 / .ZD 


1 fi1 OR QQ 
+ 01 ZD OO 


n 

y 


in to (c\ r\A\ 
lU.ly 


12.61 


(0.07) 


11.68 


(0.09) 


13 


6 


(0.3) 


ncld star 


J23134708+6127090 


OQ 


1 Q A *7 OS 
lo 4 / .Uo 


_L_fi1 07 no 

-hoi zt uy 


u 


D.vJi [y.yJo) 


16.57 


(0.17) 


20.60 


(0.16) 


29 


1 


(0.4) 




J23134689- 


-6126435 


23 


13 46.89 


+61 26 43 


5 


0.59 (0.01) 


2.24 


(0.02) 


3.11 


(0.05) 


3 


5 


(0.3) 


field star 


J23134622- 


-6126327 


23 


13 46.22 


+61 26 32 


7 


0.81 (0.01) 


1.34 


(0.02) 


1.57 


(0.04) 


3 


7 


(0.3) 




J23134589- 


^6126437 


23 


13 45.89 


+61 26 43 


7 


1.55 (0.01) 


4.99 


(0.04) 


7.45 


(0.06) 


7 





(0.3) 


field star 


J23134498- 


r-6126498 


23 


13 44.98 


+61 26 49 


8 


<0.05 


5.06 


(0.03) 


32.23 


(0.15) 


62 


8 


(0.6) 


NGC 7538 S 


J23134459- 


-6127148 


23 


13 44.59 


+61 27 14 


8 


48.09 (0.32) 


55.39 


(0.27) 


102.8 


(0.6 ) 


231 


7 


(1.1) 


nebulous object 


J23134400- 


-6126578 


23 


13 44.00 


+61 26 57 


8 


30.67 (0.25) 


28.32 


(0.14) 


170.8 


(0.7 ) 


541 


2 


(2.1) 


IRS 11 


J23134393- 


-6126503 


23 


13 43.93 


+61 26 50 


3 


3.05 (0.02) 


27.64 


(0.13) 


113.5 


(0.4 ) 


237 


8 


(1.1) 


IRS US, H 2 maser 


J23134363- 


-6126086 


23 


13 43.63 


+61 27 08 


6 


1.73 (0.04) 


5.94 


(0.04) 


8.10 


(0.09) 


25 


9 


(0.4) 




J23134302- 


-6126555 


23 


13 43.02 


+61 26 55 


5 


14.70 (0.08) 


14.07 


(0.07) 


109.9 


(0.5 ) 


401 


6 


(1.6) 




J23134293- 


-6127035 


23 


13 42.93 


+61 27 03 


5 


1.03 (0.01) 


4.84 


(0.03) 


14.4 


(0.1 ) 


37 


8 


(0.4) 




J23134160+6126427 


23 


13 41.60 


+61 26 42 


7 


1.91 (0.01) 


2.06 


(0.02) 


2.73 


(0.04) 


3 


9 


(0.3) 




J23134148+6126495 


23 


13 41.48 


+61 26 49 


5 


4.61 (0.03) 


5.41 


(0.04) 


6.77 


(0.07) 


8 


5 


(0.3) 


field star 



